Optical Lattices with Micromechanical Mirrors 
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We investigate a setup where a cloud of atoms is trapped in an optical lattice potential of a standing wave 
laser field which is created by retro-reflection on a micro-membrane. The membrane vibrations itself realize 
a quantum mechanical degree of freedom. We show that the center of mass mode of atoms can be coupled to 
the vibrational mode of the membrane in free space, and predict a significant sympathetic cooling effect of the 
membrane when atoms are laser cooled. The controllability of the dissipation rate of the atomic motion gives 
a considerable advantage over typical optomechanical systems enclosed in optical cavities, in that it allows 
a segregation between the cooling and coherent dynamics regimes. The membrane can thereby be kept in a 
cryogenic environment, and the atoms at a distance in a vacuum chamber. 

PACS numbers: 



Optical lattices for ultracold atoms Q are generated by 
standing laser light waves, obtained e.g. by retroreflecting a 
running wave laser beam from a mirror. The intensity vari- 
ation in the standing wave gives rise to a position dependent 
AC Stark shift, thus providing a periodic array of optical mi- 
crotraps for the atomic motion. In the description of atomic 
dynamics in far-off-resonant lattices, e.g. in deriving the Hub- 
bard dynamics of bosonic or fermionic atoms (T), the optical 
potentials are typically represented as c-number fields. The 
backaction of the atomic motion on the laser fields and spon- 
taneous emission into vacuum modes appear as weak decoher- 
ence mechanisms, in particular heating of atoms. Thermal and 
other vibrations of the mirror reflecting the laser light provide 
an additional source of dissipation. In light of recent exper- 
imental progress in the field of optomechanics 0, a new 
generation of optical lattice experiments is in sight, where the 
standing light waves are generated by reflection from a mi- 
cromechanical mirror or reflecting dielectric membrane (H, 
whose motion must be described quantum mechanically at 
low temperature, and which is coupled via the moving lattice 
to a distant ensemble of cold atoms. In a lattice with mod- 
erate detuning, the back-action of atomic motion on the light 
fields can be substantial [5]. This results in a sizable cou- 
pling of atomic motion and mirror vibrations, mediated by the 
quantum fluctuations of the lattice laser light. Below we de- 
velop a fully quantum theoretical description of such a setup 
(cf. Fig. [T]). In particular, we derive a Quantum Stochastic 
Schrodinger Equation (QSSE) and thus a Master Equation for 
a setup where a coherent coupling of the mirror to the col- 
lective center-of-mass motion of atoms in the lattice emerges, 
and where the mirror can be sympathetically cooled to the mo- 
tional ground state by laser cooling the atoms. We envision 
that the proposed system presents only one instance where the 
well-developed tools for laser cooling of atoms are used to 
provide a mechanism for sympathetic cooling of polarizable 
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FIG. 1 : A laser field impinging from the right is partially reflected 
off a dielectric membrane and forms a standing wave optical potential 
for an atomic ensemble. Vibrations of the membrane's fundamental 
mode will shift the standing wave field, shaking atoms in the opti- 
cal lattice. Conversely, oscillations of the atomic cloud (center of 
mass motion) will change the intensity of left/right propagating field 
components, thus shaking the membrane via changing the radiation 
pressure on it. The membrane can be kept in a cryogenic environ- 
ment, and the atoms at a distance in a vacuum chamber. 

particles (here the dielectric membrane). 

Model: Consider the setup shown in Fig.[T^. A laser (fre- 
quency lji) is shined on and partially retro-reflected from a 
dielectric membrane to form a standing wave that acts as a 
one-dimensional optical lattice trap for a cloud of atoms. The 
Hamiltonian for the system is H^u = H + Hsd, where H 
contains in a one-dimensional treatment the coupling of atoms 
and membrane vibrations to the paraxial field along the z di- 
rection as shown in Fig. [I] 
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while coupling to all other field modes is subsumed in H 3 d- 
The first two terms in Eq. ([T]) represent the energy of the fun- 
damental vibrational mode of the membrane with frequency 
uj m and annihilation operator a m , and the kinetic energy of 
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FIG. 2: (a) R- and L-modes for fields impinging from the right and 
left, respectively. All but the laser driven (cjj, i?)-mode are initially 
in vacuum, (b) Mediated action of atomic COM fluctuations on the 
membrane: At the advanced time t + atoms absorb a photon from 
the right propagating component of the laser field (with amplitude 
yfta, straight arrow) and reemit a sideband photon to the i?-field 
(fr^(t + ), wavy arrow), receiving a momentum kick of 2hki, cf. the 
term in Eq. ( [2a] ). At a later time t the sideband photon is annihilated 
0>R(t)) and transferred back to the laser field giving a momentum 
kick to the membrane, cf. first term in Eq. (2b) . (c) Reverse action: 
Via corresponding emission and reabsorption processes as given in 
Eqs. pb| ) and j2c] > membrane vibrations affect atomic COM fluctu- 
ations. In these processes both the R- (upper panel) and the L-field 
(lower panel) contribute, effectively mediating a stronger action than 
the process shown in (b). 

atoms with momentum pj and mass m at . The third term is 
the level shift operator [6 ] describing the dispersive interac- 
tion of the N two level atoms (with dipole matrix element 
H) with off-resonant field modes with detuning 5 from the 
atomic resonance within a frequency band $ <C 5 Q. In 
this treatment electronically excited states are already adia- 
batically eliminated. The second line in Eq. ([]]) is the po- 
tential of the radiation pressure force on a slab of thickness 
I, cross section A, index of refraction n, and a small excur- 
sion Z of the membrane from the equilibrium. The radia- 
tion pressure force is proportional to the difference of the in- 
tensities on either side (at ±|) of the membrane. The elec- 
tric field is E(z) = E R (z) + El{z) and consists of the 
two ID continua of plane wave modes impinging from the 
right (R) and left (L), respectively, cf. Fig. [2^. The ex- 
plicit mode functions can be found in 0, and will be de- 
noted by A a (k,z) (a = L, R) consistent with the notation 
used there. Thus, E+(z,t) = £ JdujA a (k, z)b U:0i e~ iujt with 
[fru;,a> fyji> ^] = S a ^5(uj — w) and k = uj/c. We are using 

a narrow band approximation £ = \J a ^ € 1 qA and take the 
fields in an interaction picture with respect to the free field 
Hamiltonian. The contribution of other than paraxial modes 
contained in H 3 d can be largely suppressed in the follow- 
ing. Its effect will be relevant and addressed separately in the 
context of momentum diffusion of atoms in the lattice. The 
driving laser of frequency uji is at this point included by as- 
suming the field to be in a coherent state \a) = D|vac) (with 
D = exp(a*b L j li R — h.c.)), corresponding to a laser field of 
power P = huia 2 /27r (with photon flux a 2 /2ir) impinging 
from the right. 

Quantum Stochastic Schrodinger Equation: We are inter- 
ested in treating those processes only whose amplitudes are 
enhanced by the laser amplitude a in either first or second 



order. In order to identify those we apply a Mollow transfor- 
mation ifTOl and move to a displaced picture, where the field 
is effectively in vacuum |vac) = D^\a) and the new Hamil- 
tonian is accordingly H' = D^HD. In H' the second or- 
der terms in a provide the lattice potential (with trap depth 

Vo = 4/xWVi an d trap frequency u at = for a 

membrane power reflectivity r), and an unimportant mean ra- 
diation pressure force on the membrane. The terms of first 
order in a, which contain the coupling of the motion of mem- 
brane and atoms to electromagnetic (EM) field modes ini- 
tially in vacuum, are treated in a Lamb-Dicke expansion in 

(h&at) where i at — m ^ - . The result can be interpreted 
as a (Stratonovich) quantum stochastic Schrodinger equation 
(QSSE) (7) with time delays 

£\V) = -iH(t,t-,t+)\V) 

= | - iH sys + x at (^/g at , R b R (t + ) - h.c.) (2a) 

+ (iy/9^Rb R (t) - ^gm,Lb L (t) - fl.C.) (2b) 

- x a t{y/9^Rb R (t-) + iy/g^Eb L {t~) - h.c.)}|*), (2c) 

where H sys = (jJm(4n a m + u at a} at a at denotes the sys- 
tem Hamiltonian for the free evolution of the membrane 
and the atomic center of mass (COM) mode, with annihila- 
tion operator a at for the later. The remaining terms denote 
the coupling of position fluctuations £ m (at) = ( a m(at) + 
h.c.)/\/2 of these modes to vacuum fluctuations of the EM 
field described by slowly varying field operators b a (t) = 
J -^Lb a {uj) exp[— i{uj — uji)t] (a = i?, L), evaluated at time 

t as well as at advanced and retarded times ^ = t ± ^, 
where d is the distance between the atomic ensemble and the 
mirror, cf. Fig. [2]). In this Hamiltonian the linear extension 
of the ensemble is neglected and we used the approximation 
b(t ± z -f) ~ b(t ± f ) for all atoms fl]Q. The coupling con- 
stants of membrane vibrations and atomic COM position fluc- 
tuations to the EM field are given by 

g m , R = 2(ak l i m ) 2 x 2 jix, g at , R = 27rN(cj at /4ak l £ at ) 2 , 

where £ m = y m h u (rn rn is the effective mass of the mem- 
brane), and g m (at),L = | 9m(at),R with the membrane trans- 
mittivity t = 1 — r. In Fig.[2j) and c we provide an intuitive in- 
terpretation of each term in the QSSE ^ as an (Anti)Stokes- 
backscattering process, where a laser photon is absorbed or 
emitted together with, respectively, an emission or absorption 
of a sideband photon at frequency uoi ±u m ^ at ^ in one of the ini- 
tially empty modes, along with a momentum transfer of 2hk\ 
to either the membrane or the atomic COM mode. The effec- 
tive mediated coupling between atoms and membrane, as well 
as radiation pressure noise on each system, will eventually re- 
sult from second order processes, where sideband photons are 
emitted and reabsorbed, Fig. (2]} and c. 

Cascaded System Master Equation The corresponding ef- 
fects of these second order processes can be seen explicitly 
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when Eq. ^ is transformed into an equivalent master equa- 
tion (MEQ) without time delays in the limit —> t, and under 
the Born-Markov approximation for the EM field in vacuum. 
The technicalities of this procedure, which requires some care 
regarding correct time ordering, are outlined below. The result 
is a master equation for the density operator p of the mem- 
brane and the COM motion of the form 

P = -i[Hsys + gXatXm, p] + Cp + L m p + L at p. (3) 

This MEQ contains, firstly, an effective atom-membrane inter- 
action with a coupling g = 2(^g rn , R g at)R + y/g m ,L9at,L) = 

uj at ' N ™ at • Here we assume comparable oscillation frequen- 
cies uj m uj at . While the ratio of masses is always 

small under reasonable conditions, the rate of coherent cou- 
pling g can still be significant thanks to the collective enhance- 
ment by y/N. 

Secondly, we find a contribution 

Cp= ^Tjr ([%rm x atp] ~ [p x at 7 x : ra]) •> 

which is proportional to the membrane transmittivity t. The 
peculiar form of this term is in fact generic for cascaded 
quantum systems Q HS1 - Its effect becomes clear e.g. when 
looking at the evolution of mean values according to Eq. ([3]), 

(fiat) = g(xm)^ while (pm) = gt(x at ) H , where dots 

stand for contributions of H sys ,L m and L at . Thus, the action 
of atoms on the membrane is reduced by a factor t as com- 
pared to the action of membrane fluctuations on atoms lfT9l . 
This can be understood from the fact that action of membrane 
on atoms is mediated through both fields, R and L, while only 
the R field contributes to the action of atoms on the mem- 
brane, see Fig. [2] 

In the limiting case t = 1, that is for an ideal mirror, 
the term Cp does not contribute and the interaction is only 
due to the Hamiltonian part of Eq. ([3]). The scaling of g 
as given above can then be understood in terms of a sim- 
ple consideration in terms of quasi-static field modes: The 
mirror displacement Z(t) provides a time-dependent bound- 
ary condition for the EM field, such that the standing wave 
laser field is E(z) oc sin[fc/(z — Z)). The corresponding lat- 
tice potential for the atoms is V{zi) = ^ Vb sin[fc/(^ — 
z )? - T,iVotf(Szf + 25zi5Z), where 5z 4 and 5Z de- 
note small fluctuations of atoms and mirror around their re- 
spective equilibrium positions. This is on the one hand just 
the harmonic potential for the atoms with a trap frequency 
uo at = (2Vokf /mat) 1 / 2 (the corresponding term for the mir- 
ror can be absorbed into cj m ), and on the other hand a direct 
coupling of the mirror and the COM motion of the atoms with 
a strength 2Vokf J2i SziSZ — hg x at x m, correctly reproduc- 
ing g = VN 2V k?£ at £ m /h = uj atx J^^. However, for the 
realistic case of a membrane with reflectivity t < 1 such con- 
siderations can hardly lead to the correct master equation ^ 
(2(2. 

The strength of coherent coupling has to be compared to 
various rates of decoherence, described in Eq. ^ by a num- 



ber of Lindblad terms L m , L at . One channel of decoher- 
ence is caused by coupling to the EM vacuum field, which 
is contained in the QSSE in Eq. ^ and results in radia- 
tion pressure induced momentum diffusion described by Lind- 
blad terms L d $ m) p = htl{L) D l x at(rn)]p, where D[a]p = 
2apo) — a^ap — pa^a. The atomic momentum diffusion rate 
i s 7a = (Wat) 2 7se^b/^ with a natural linewidth j se and 
a detuning S from atomic resonance. Note that a momentum 
diffusion rate j^t^ f° r me atomic COM motion is equal to 
the well-known single atom diffusion rate 1 8 ] only if all inter- 
actions between atoms are neglected. In the ID treatment of 
the QSSE ^ these interactions are in fact overestimated, as 
they do not fall off with the distance between atoms. Only a 
more careful 3D treatment, taking into account all terms ne- 
glected in Hsd, reveals the correct distance dependence of 
dipole-dipole interactions and allows to ignore these effects 
for a dilute sample. The diffusion rate of the membrane mode 
due to radiation pressure is jfflf = (4xP / m m c 2 )(uoi / uj m ) 
(where P is the laser power). This last result is consistent 
with calculations in lfT3l on the effects of radiation pressure 
on a mirror in free space. 

On top of these decoherence channels there will be ther- 
mal heating of the membrane due to clamping losses and 
laser absorption, which drives the vibrational mode under 
consideration to thermal equilibrium at temperature T, i.e. 
L mP =1 f(n + l)D[a m ]p + ^nDlal^p, with thermal oc- 
cupation ft — kBT/huj m , and 7 m = uJm/Qm for a mechan- 
ical quality factor Q m . For atoms in turn, it is important to 
note that the dissipation processes L at can to a large extent 
be engineered. In particular, it is possible to use Raman side- 
band laser cooling to cool atoms individually close to their 
ground states in the optical lattice potential ifTTIl . In the mas- 
ter equation we account for this by adding a Lindblad term 
for the COM mode, L c a ° t ol p = \^ c a °° l D[a at ]p with a Raman- 
cooling rate . If laser cooling is switched off, only the 
background decoherence at a rate remains. The dynam- 
ics of the master equation ^ thus provides two interesting 
regimes: with efficient laser cooling of atoms it is possible to 
sympathetically cool the mirror, while for <C g a regime 
of coherent evolution becomes accessible, limited only by the 
decoherence rate of the membrane 7^ = 7 m n. 

Case Study We estimate the magnitude of the relevant 
processes as follows: For a small SiN membrane of dimen- 
sions 150 pm x 150 pm x 50 nm exhibiting a mechanical res- 
onance at uj m = 2tt x 0.86 MHz with an effective mass 
of m m = 8 x 10 _13 kg and a high mechanical quality of 
Q m = 10 7 at T < 2 K, a power reflectivity t = 0.31 can be 
achieved for the wavelength of A = 780 nm relevant for 87 Rb 
l4l . A lattice beam with power P = 7mW and a waist of 
230 pm at the position of atoms, detuned by S = 2tt x 1 GHz, 
can provide a sufficient lattice depth such that the longitudinal 
trap frequency uo at ~ uo m . For an ensemble of N = 3 x 10 8 
atoms IfTTIl a strength of coherent coupling g = 40 kHz can 
be achieved. For the decoherence due to radiation pressure 
noise one finds a momentum diffusion at rate 7^^ = 52 Hz 
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for the membrane, and an atomic momentum diffusion in the 
lattice at rate = 16 kHz. For the membrane the dom- 

inating dissipative effect will clearly be thermal decoherence 
at rate 7^ = 4 MHz (24 kHz) at room temperature (2 K). On 
the other hand, Raman sideband cooling of atoms at a fast 
rate ^ ol ~ 20kHz is possible fiTTl . For these parameters, a 
regime where u m = u at > g > -f c a ° t ° l ~ 7^ > 7^ /J > 
7^ is accessible. To suppress trap loss due to light-assisted 
collisions, blue detuning of the lattice laser is advantageous 
|[TT1l . Transversally, the atoms can be confined in a far-detuned 
2D lattice, so that atomic densities of order 10 12 atoms/cm 3 
are realistic lfT2ll . A concern is the spread of vibrational fre- 
quencies Auj at ~ 2tt x 24 kHz across the ensemble due to 
the Gaussian intensity profile of the laser. It can be reduced 
by using a beam with a flat top profile. Due to its modularity, 
the presented setup is very flexible, e.g. several atomic ensem- 
bles could be trapped in multiple foci of the same lattice laser, 
thereby enhancing N and thus g. 

Sympathetic Membrane Cooling In order to evaluate the 
corresponding efficiency of sympathetic cooling of the mem- 
brane via laser cooling of atoms, we solve the master equa- 
tion ([3]) for the steady state occupation n ss = (aj n a m ) ss . For 
the parameters given above one finds a cooling factor / = 
n/n ss ~ 2 x 10 4 , which is sufficient for ground state cooling 
(n ss ~ 0.8) of the membrane when starting from 500 mK. In 
order to get some more insight into the cooling efficiency we 
consider the weak-coupling limit ^> g. In this case we 
can eliminate the atomic COM mode adiabatically along the 
lines of the treatments of the equivalent problem of optome- 
chanical laser cooling [21 13. For the mean occupation one 



finds a rate equation ^ (aj^ 



an effective cooling rate 



7m + 



-r m ((a] n a m ) -n ss ) with 
and a final occupa- 



tion n ss ^^71+ ("2jjj— ) • For large enough cooling rate the 
thermal contribution can be suppressed and the limitation is 
only due to Stokes-scattering processes. As in laser cooling of 
ions or optomechanics, this can be suppressed in the resolved 
sideband limit, which amounts here to having <C uj m . 

Under this condition ground state cooling is possible. 



If laser cooling of atoms is switched off (7; 



cool 



0), a 



regime of coherent coupling is accessible, at least for cryo- 
genic temperatures, where g > 7^,7^^- As compared to 
the usual optomechanical setup, where the equivalent param- 
eter to Jat° l * s me cay ity decay rate which is a fixed parameter, 
this is a qualitatively new feature of the setup considered here. 
It is well known, and has been extensively studied in various 
contexts ifTTlL that the given coupling Hamiltonian ~ x at x m - 
which becomes ~ (a m a^ t + a^a a t) in the rotating wave ap- 
proximation {uj m ^> g) - allows for a coherent state exchange 
of the two systems. 

Appendix: Conversion of the QSSEfrom Stratonovich to Ito 
form We still need to explain how to transform the QSSE 
with time delays to the Markovian master equation ([5J. We 
first transform the Stratonovich QSSE to an Ito QSSE by first 
integrating Eq. ([2J up to a time t + At (such that At ^> r ^> 
l/fi, where $ is the bandwidth of modes) and then expand to 
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FIG. 3: Cooling factor / = n/n ss versus effective coupling g be- 
tween the atomic COM motion and the membrane vibrations and 
rate of Raman sideband laser cooling 7a? oZ , as determined from an 
exact numerical solution of the effective master equation {3} for a 
resonant system cj m = uo a t and a mode of the membrane with 
uj m = 2tt x 0.86 MHz and Q m = 10 7 at T = 2K. The calcula- 
tion includes radiation pressure induced momentum diffusion at rate 



7, 



diff 



16 kHz for atoms and 7™ 7 = 52 Hz for the membrane. 
The parameters g ~ 40 kHz and jat° l — 20 kHz discussed in the 
text provide a cooling factor of 10 4 , and yield a steady state occupa- 
tion n ss below one if the system is precooled to below 500 mK. At 
room temperature the sympathetic cooling effect will be still clearly 
observable. 



second order, 

tf(At)|*(t)> = [l-if^dhHitut^tt) (4) 
" It +At dt i It" dt 2 H(U, tl, t+)H(t 2 , t~ 2 , t+)} |*(t)> 

R9at,R + ^9m y L9at y L)XmX a t At 

- l(gm,R + 9m,L)x 2 m At- \g at , L x 2 at At}\^ (t)) . (5) 



We use here the Ito increments AB a (t) = f* +At dt'b a (t') 
and their property AB a (t) \ \£) = for all times. In the double 
time integral care needs to be taken in dealing with the time 
delays. It is easy to check that J^ +At dt\ f* 1 d^fra (£1)^(^2+ 
S) = 9 (5)5 a , p At I2TII . After applying this rule it is possible to 
take the delay to zero, t^ — » t, as has been done in the second 
step of Eq. ([?]). The last equation ([5]) can then be interpreted 
as an Ito QSSE, which can in turn be converted to the mas- 
ter equation ^ following standard procedures |7|. Note that 
this procedure results in a seemingly collectively enhanced 
atomic momentum diffusion, which is an artefact of the ID 
model adopted in Eq. ([T]). A more careful treatment based on 
a 3D model for the coupling of atoms to the EM field cor- 
rectly yields the well known individual momentum diffusion 
of atoms in an optical lattice at the rate given in Eq. ([3]), along 
with dipolar interactions which are suppressed at low densities 
If22l 
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